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We extend an earlier semiclassical model to describe the 
dissipative motion of atoms coupled to M modes inside a 
coherently driven high-finesse cavity. The description includes 
momentum diffusion via spontaneous emission and cavity de- 
cay. Simple analytical formulas for the steady-state temper- 
ature and the cooling time for a single atom are derived and 
show surprisingly good agreement with direct stochastic sim- 
ulations of the semiclassical equations for A'^ atoms with prop- 
erly scaled parameters. A thorough comparison with standard 
free-space Doppler cooling is performed and yields a lower 
temperature and a cooling time enhancement by a factor of 
M times the square of the ratio of the atom-field coupling 
constant to the cavity decay rate. Finally it is shown that 
laser cooling with negligible spontaneous emission should in- 
deed be possible, especially for relatively light particles in a 
strongly coupled field configuration. 

PACS numbers: 32.80. Pj, 33.80.Ps, 42.50.Vk 



I. INTRODUCTION 

For a large variety of atoms, laser cooling has proven to 
be an extremely successful and widely applicable method 
to generate fairly large ensembles of ultra-cold neutral 
atoms 1^. Nevertheless, some intrinsic effects, like re- 
absorption of the spontaneously scattered light, optical 
pumping to unwanted states or dipole heating, limit its 
applicability and the final phase space density that can 
be achieved. It has been suggested that these prob- 
lems could be partly cured by the help of high-finesse op- 
tical resonators to enhance the cooling process and make 
it more selective in frequency space as well as in posi- 
tion space. For single atoms moving in the field of only 
one or a few photons in high-Q microscopic resonators 
these results could largely be demonstrated experimen- 
tally 1^,0. Although this nicely confirms the theoreti- 
cal concepts and model approximations, the scaling to- 
wards larger particle and photon numbers as well as large 
interaction volumes is not entirely clear. Some initial 
results on the scaling properties of this cooling scheme 
were obtained for ensembles of well trapped particles [D , 
where a harmonic approximation for the trapping poten- 
tial can be assumed. Neglecting cavity induced momen- 
tum diffusion, collective equations for the atomic motion 
could be derived from which the scaling properties with 
atom number and volume could be read off. Similarly in 
the very far detuned weak field case, where spontaneous 
emission and field induced atom-atom interaction can be 
neglected, estimates of the scaling properties could be 



derived from geometrical considerations ||^. 

The range of the validity of both of these approaches 
is, however, not entirely clear and does not cover a large 
parameter regime. In many practical cases neither the 
tight binding situation allowing a collective description 
of the properties of a whole ensemble, nor the weak field 
assumption treating each atom separately will be real- 
ized. In that case one cannot expect reliable answers from 
these models. Only a more general description bridging 
the gap from the single atom to the collective descrip- 
tion and from a single photon to a coherent field with 
many photons should provide adequate answers on the 
efficiency and realizability of this cooling scheme. 

In this work we use a quite general semiclassical ap- 
proach which only relies on the assumptions of a weak 
saturation and a not too low temperature T 3> TrecoU, 
which should be true for most practical setups. The 
model as defined in section II consistently includes spon- 
taneous emission, cavity decay, atom-atom interaction 
via the field, and realistic mode functions. For a single 
atom in a standing wave field simple analytic expressions 
for the final temperature and the cooling time are de- 
rived in section III. In a rescaled version these are then 
applied to the A'^-atom case and tested against numerical 
simulations in section IV. In section V we compare the 
results to standard Doppler cooling and discuss the possi- 
bility of cooling without spontaneous emission in section 
VI. Finally, in section VII we briefly discuss the scal- 
ing properties for the case of several degenerate cavity 
modes. 



II. SEMICLASSICAL MODEL 

We consider a dilute gas of N identical two-level atoms 
of mass TO, transition frequency uja, and spontaneous 
emission rate F interacting with M modes of a high- 
finesse optical cavity with resonance frequencies ujk , cav- 
ity decay rates Kk, and mode functions fk{x), k — 1..M, 
respectively. The mode functions are mutually orthogo- 
nal and fulfill the normalization condition j \ fk\'^ dx = V , 
where V is the mode volume. The spatially averaged sin- 
gle photon Rabi frequency is denoted by g. The cavity is 
driven by a pump laser of frequency LOp, the spatial over- 
lap of the pump with the mode functions at the driven 
cavity mirror yields pump strengths rjk. 

Generalizing previous work we derive a semiclas- 
sical model of the system described by a set of coupled 
stochastic differential equations (SDEs) for the atomic 
positions a;„, momenta p„, and mode amplitudes afc. 
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where 5(a;) = X^fe fk{x)ak is the field ampHtude at po- 
sition X, [/q = Aa5'^/(A^ -|-r^) the hght shift per photon 
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j^)^ 7 = r,g^/(A^+r^) the photon scattering 



rate, and A^: ~ Up — Uk the detuning of the /cth mode 
from the pump laser. 

The interpretation of the various terms in eqs. (|l]) is 
rather simple. The two terms in eq. ( p^ ) correspond to 
the dipole force and the radiation pressure force, respec- 
tively. Equation (|l^) describes the pumping of the cav- 
ity modes, the frequency detuning from the pump laser 
shifted according to the coupling to the atoms, and the 
mode damping due to cavity decay and photon scatter- 
ing by the atoms. dP„ and dAk are white noise incre- 
ments which can be characterized by a diffusion matrix 
Dij = (dFidFj), where Fi = {Pn,Ak}. The diffusion 
matrix exhibits nontrivial cross correlations between mo- 
mentum diffusion and cavity amplitude and phase fluctu- 
ations pt] . We may decompose the noise terms by writing 



dP„ = dP„*f + dP, 
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The first terms correspond to momentum and mode am- 
plitude fluctuations due to spontaneous emission of pho- 
tons into vacuum modes. The second terms describe the 
correlated momentum and field fluctuations, which oc- 
cur, when the nth atom scatters a photon among differ- 
ent cavity modes. The precise form of all noise terms is 
given in Appendix ffl^ 

The set of SDEs (^ is well suited for a numerical im- 
plementation. Note that the number of equations grows 
linearly with the number of atoms and modes, in sharp 
contrast to full quantum descriptions, which become ex- 
ceedingly large for all but the simplest situations. The 
inclusion of the quantum noise terms, on the other hand, 
allows to obtain reliable results, for example, for steady- 
state temperatures and trapping times, which cannot be 
derived from a simple classical model. 

Let us now specify the parameter regime we are con- 
centrating on in this paper. We assume an optical cavity 
of sufficiently high finesse with large separation between 
the various eigenfrequencies. Then we can neglect all 
modes, which are not close to resonance with the pump 



field frequency. For a standing wave cavity the param- 
eters can be chosen such, that only a single mode will 
fulfill this condition. For a ring cavity setup one should 
include at least the two degenerate counterpropagating 
running waves, while in a quasiconfocal or quasispherical 
setup many transverse modes can be involved in the dy- 
namics. We will discuss these situations later in section 
VII. For the following, we will specialize the model to 
a one-dimensional standing wave setup, neglecting trans- 
verse atomic motion. Therefore, only a single mode func- 
tion f(x) — cos(fcx) is considered and we will drop the 
subindex k in the notation for the field amplitude, de- 
tuning, and pump strength. Since the derivation of the 
SDE is based on an adiabatic elimination of the internal 
atomic degrees of freedom, we must ensure a small atomic 
saturation parameter, that is, s = |apg^/(A^ -l-F^) ^ 1. 
To this end we assume a large atomic detuning Aq 3> F. 
In order to observe a significant effect of the atomic mo- 
tion onto the cavity dynamics, we require A « —n and 
-^|C^| ^ 1^- Hence, for large atom numbers N, we assume 
|J7o| 1^- Although these parameters are known to give 
rise to cooling independent of the sign of the atomic de- 
tuning (and hence the sign of Uq), we will in the following 
use negative values of Uq since, in a full model, this ad- 
ditionally gives rise to three dimensional trapping. 



III. ANALYTICAL RESULTS FOR A SINGLE 
ATOM 

Before we discuss the results obtained from numerical 
simulations of the SDEs, we will briefly summarize the 
analytical results obtained earlier for a single atom in 
a single mode of a standing wave cavity and specialize 
these to the present parameters. As we will see later, 
this provides an excellent basis for the discussion of the 
multi-atom case. 

The analytical results [§,0 were derived under the as- 
sumption of weak pumping r] such that there is at most 
one excitation in the system. However, it can be shown 
that under the conditions required for the adiabatic elim- 
ination of the atomic excited state, these results also hold 
for higher photon numbers as long as s ^ 1. Thus, we 
may apply the formulas in the semiclassical limit dis- 
cussed here. For atomic velocities v fulfilling kv < k the 
force on a moving atom can be expanded to lowest order 
in kv/n as 



F = Fo + Fiv, 



(3) 



where Fq is the stationary dipole force from eq. ( p^ for 
a fixed atomic position, and Fi is the linear friction co- 
efficient, 

Fi{x)^~2k^cos{kxfsm{kxf'^^. (4) 

Neglecting localization effects, the position-averaged fric- 
tion force thus reads 
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Similarly, the averaged momentum diffusion coefficient is 
found to be 



D = k^K 
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and hence the steady-state temperature is 
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Thus, we can get sub-Doppler cooling for k <V . 

Finally, from the friction coefffcient we obtain an esti- 
mate for the cooling time defined from the exponential 
decrease of the kinetic energy E{t) of the form 



E{t) 



exp(-t/Tc) -I- 
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with 



2|Fi| v^U^ 
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Here ujji — fik"^ /{2m) is the recoil frequency. Hence 
we get faster cooling for larger cavity field intensity and 
larger optical potentials per photon. Moreover, the cool- 
ing time is proportional to the atomic mass such that the 
scheme works best for relatively light atoms. 



IV. COOLING OF MANY ATOMS 

Let us now consider N atoms. As they all couple to the 
same mode, one might expect that in general they will 
perturb each other and interact via this mode [[ll|--[l4[. 
In fact, for a few atoms and a few photons inside a 
high-finesse cavity such effects have recently been demon- 
strated experimentally [ p^ . In order to investigate the 
cooling properties of many atoms, we numerically inte- 
grate the SDEs and average over many different real- 
izations. Choosing a flat initial distribution of atoms in 
space and a large initial kinetic energy we calculate the 
time evolution of the kinetic energy. In the following we 
will concentrate on the steady-state temperature and, in 
particular, on the cooling time scale Tc. 

In Fig. y we plot the numerically obtained mean ki- 
netic energy as a function of time for a single atom and 
for ten atoms. For the many atom case we rescale the pa- 
rameters such that NU{) and the optical potential depth 
(proportional Uqi^ / k?) are constant. We find that the 
curves are very well approximated by an exponential fit 
as in eq. (||) with cooling times Tc = 142n~^ {N = 1) and 
Tc = IIIOk"^ {N = 10). For the given parameters, the 
analytical estimate given by eq. (||) yields Tc = 118.6k^^ 
and 1186/t~^, respectively, which is surprisingly good. 
Therefore, the single atom cooling time appears to make 



sense for the many atom case. However, the operating 
conditions differ for varying atom numbers, according to 
the condition that NUo remains constant. 
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FIG. 1. Time evolution of the mean kinetic energy per 
atom obtained from numerical simulations of the SDEs aver- 
aged over 100 realizations (solid curves). The dashed curves 
are exponential fits. Lower curves: A'^ = 1, upper curves; 
N = 10. The parameters are NUq = -0.6k, N'^j = Q.03k, 
A = -Q.6k, 7? = sVNh, k = 4:15ujR. 
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FIG. 2. Steady-state kinetic energies per atom (a) and 
cooling times Tc (b) for atom numbers from 1 to 100. Crosses 
show the numerical results for A — —0.6k, circles for A — —k. 
All other parameters are the same as in Fig. ^ The solid and 
dashed lines are linear fits to the numerical data. 
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The steady-state kinetic energies obtained from the ex- 
ponential fits of Fig. |l| are 466^^? {N = 1) and SlOw/j 
{N = 10). While these are approximately of the same 
size, they are significantly larger than the value of k/4 « 
lOOcj^f obtained from eq. However, cq. (j^) was de- 
rived assuming Uq ^ k and a fiat spatial distribution in 
the steady state. Both of these conditions are violated 
for the present parameters. 

It should also be emphasized that an exponential fit 
as in Fig. Q only works for sufficiently hot atoms where 
the friction force (^) can be replaced by its mean value 
and can therefore be treated as constant. For very cold 
and well localized atoms the position dependence of the 
friction must be taken into account. In fact, eq. (||) pre- 
dicts that the friction force vanishes completely at the 
bottom of the potential wells. Hence, if in the course of 
the cooling process the localization increases, the time 
scale of the cooling increases and the kinetic energy no 
longer decreases according to an exponential law. For 
such parameters, a unique cooling time Tc can only be 
deduced from the initial stages of the cooling while the 
atoms are still uniformly distributed in space. 

For the same scaling of the parameters with N as dis- 
cussed above, we plot the steady-state kinetic energy and 
the cooling time for particle numbers from one to 100 in 
Fig. 1^. We find an approximately constant final temper- 
ature and a cooling time proportional to N. Taking the 
scaling of Uq with the particle number N into account, 
these results agree with the predictions of the single-atom 
estimates of the previous section. 

In Fig. ^, a different scaling of the system parameters 
is applied. We choose a small value of Uq and keep all pa- 
rameters constant while changing the atom number. For 
the steady-state temperature we obtain approximately 
constant values as long as iV|J7o| ^ k/2 {N ^ 10). For 
10 ^ < 25 the frequency shift induced by the atoms 
is of the order of the cavity detuning A = — k, which 
reduces the cooling efficiency and consequently increases 
the steady-state temperature. For even larger atom num- 
bers the cavity is shifted into positive detuning where the 
friction force changes sign and the atoms are accelerated. 
In this limit no stationary state is achieved and no cool- 
ing time can be defined. 

For the same parameters, figure ^(b) shows the cooling 
time for varying atom number. We see that Tc slightly 
decreases with increasing N. Close inspection of the nu- 
merical data shows, however, that this decrease can be 
entirely explained by an increase of the cavity photon 
number according to the increase of NUq, which shifts 
the cavity closer to resonance for larger atom numbers. 
Keeping the photon number constant, e.g., by adjusting 
the pump strength rj, would thus give rise to a cooling 
time essentially independent of N. 

Therefore, the numerical simulations show that for 
NUq < K the cooling time for an ensemble of A'' atoms 
is of the same order of magnitude as the cooling time for 
a single atom. This suggests that the individual atoms 
in the cloud are cooled independently from each other, 



although they are all coupled to the same cavity mode. 
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FIG. 3. Steady-state kinetic energies per atom (a) and 
cooling times Tc (b) for atom numbers from 1 to 10. Pa- 



rameters are Uo = 
K — 415cJij. 
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For the case of NUq <C k (no effect of the atoms on 
the cavity mode) this has also been predicted by Vuletic 
and Chu . The basic argument there is that the cavity 
enhances backscattering of photons by moving atoms ac- 
cording to the Doppler shift, similar to free-space Doppler 
cooling. We now see that the cooling of atoms is nearly 
independent of the atom number also in cases where the 
cavity photon number is significantly changed by the 
presence of the atoms {NUq ^ k). We attribute this 
to the fact that, considering the motion of an individual 
atom out of the whole ensemble, only the change in the 
cavity field induced by that atom itself is correlated with 
the atomic motion, while the influence of all other atoms 
is uncorrelated and hence on average cancels. Note that 
this argument is based on the assumption of no corre- 
lation between the motion of different atoms. Thus, for 
deep (harmonic) optical potentials in the classical limit 
where all particles oscillate with the same frequency a 
completely different scaling with the atom number has 
been found For the parameters considered here, the 
anharmonicity of the sinusoidal potential implies a dif- 
ferent oscillation frequency for different atoms. Together 
with the independent part of the momentum diffusion 
this ensures fast dephasing between the individual atomic 



4 



oscillations. Despite the fact that they commonly inter- 
act with the same field these mechanisms guarantee that 
the atoms can be considered as independent. Such decor- 
relation (phase space mixing) has also been shown to be 
important for stochastic cooling of atoms in optical traps 
1^,0, which has been suggested as another alternative 
to cool atoms without spontaneous emission. 



V. COMPARISON WITH FREE-SPACE 
DOPPLER COOLING 

Since we have seen that the single atom formulas of 
Sec. [I] provide a good approximation for the many atom 
case, too, we will investigate these in the following a bit 
further. 

Let us first introduce the spatially averaged atomic sat- 
uration s as 



2g^cos(fcx)2 ^ ti'^Uq 



r2 



(10) 



For the validity of the SDEs (|l|) we required s ^ 1. We 
can then rewrite the averaged friction force (^ in the 
form 



and the cooling time as 



(11) 



(12) 



Hence, for a given atom and a given saturation the max- 
imum friction force and minimum cooling time Tc is lim- 
ited by the ratio of the atom-cavity coupling g and the 
cavity decay rate k. Cooling is most efficient if this ratio 
is large, that is, for high- finesse cavities. 

For standard Fabry-Perot resonators with mirror 
curvature larger than the cavity length /, the mode waist 
scales as ^^Z* and thus the mode volume V as l^^'^. Since 
g (X 1/VV and k cx 1/Z we find (g/n)^ oc Vl. Thus, 
cavity enhanced cooling seems more favorable for longer 
cavities, which also would allow to include more atoms. 
However, this route is limited by the capture range of 
the cooling given by kv < k. Thus, for a small cavity 
decay rate only already very cold atoms are subject to 
the cavity cooling force. 

Let us now compare eqs. (^l|) and (|l2|) with the corre- 
sponding results for free-space Doppler cooling ||l^ . This 
provides optimum cooling for Aq = — F. The friction 
force is given by 



pDopp ^ _^^2^ 

and the cooling time is therefore 

Dopp ^ 1-1^ 
4s ^ 



(13) 



(14) 



We note that the only difference with the cavity-induced 
cooling is the absence of the factor {g/ kY . Hence, only 
in high-Q cavities with g > k the cooling scheme dis- 
cussed here is faster than standard Doppler cooling. The 
achievable steady-state temperatures of the order of k/2, 
on the other hand, can be significantly lower than the 
Doppler limit of r/2 even for cavities with g < k but 
at the expense of a smaller capture range as already dis- 
cussed before. 



VI. SPONTANEOUS ATOMIC DECAYS 

Cavity-induced cooling has been suggested as a possi- 
ble means to cool particles without closed optical transi- 
tions such as certain atomic species or molecules. The 
motivation behind this proposal is that the cooling 
scheme relies on induced internal transitions of the parti- 
cle rather than on spontaneous emission of photons. On 
the other hand, equation (11) shows that for large atomic 
detuning and correspondingly small saturation, which is 
required for a small photon scattering rate 7, the cooling 
time also becomes very long. The figure of merit in this 
respect is the number of photons Nph which are spon- 
taneously scattered per particle during one cooling time 
Tc- Only for Nph ~ 1 effective cooling of, for example, a 
molecule becomes feasible. 

Using the results of the previous section, we can easily 
obtain an estimate for this number from the definition 
Nph = 2TsTc yielding 



Nph = 



_ 1 F 
2ljrU^ 2 LJR 



(15) 



Note that this result is independent of the pump strength 
r] and the atom detuning , and is entirely fixed by atom 
and cavity properties. 

As in the previous section, it is interesting to compare 
this with the corresponding result for standard Doppler 
cooling given by 



N 



Dopp 
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Again, the only difference is the missing factor of {n/g)^. 
Therefore, only for g > k the number of spontaneously 
emitted photons per cooling time is reduced with re- 
spect to free-space Doppler cooling. Despite the fact that 
cavity-induced cooling works at arbitrarily large atomic 
detuning, the increased cooling time counteracts this re- 
duction of the scattering rate 7, such that the over- 
all efficiency falls behind Doppler cooling for cavities of 
not sufficiently high finesse. For high-Q cavities, on the 
other hand, the reabsorption problem, which limits the 
achievable particle density in free-space laser cooling, is 
strongly reduced too, as a photon which has decayed 
through the cavity mirror cannot be reabsorbed po|| . 
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VII. MULTI-MODE CAVITIES 

Let us now return to the more general situation where 
the optical cavity supports several (nearly) degenerate 
modes as described by the SDEs (|l|). 

In a first step, we use the friction coefficient obtained 
for a single atom in a ring cavity [2l|] to derive analytical 
estimates analogously to those in sections III and VI for 
a single-mode standing-wave cavity. We findFI 





(17) 




(18) 


1 k4 ^ 


(19) 
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(20) 



Hence, for the same mean atomic saturation we obtain 
the same steady-state temperature but only half the cool- 
ing time and only half the number of spontaneously scat- 
tered photons per atom compared to the single- mode cav- 
ity. 

This suggests that using M degenerate optical modes, 
for instance in a near-confocal or near-planar setup, im- 
proves the efficiency of cavity-enhanced cooling by a fac- 
tor of M. In fact, numerical simulations of ten atoms 
and one to ten modes have confirmed this simple scal- 
ing. Therefore, a multi-mode cavity reduces the cooling 
time and the number of incoherently scattered photons 
per atoms by a factor of M{g/ nf' compared to free-space 
Doppler cooling. Finally, it should be emphasized that 
in this case and using again A « — k, cooling occurs if 
MN\Uq\ (instead of Af|C/o| for a single mode) is of the 
order of or smaller than k. 



VIII. CONCLUSIONS 



the atom number. The analytical results obtained ear- 
lier for that situation can thus be specialized to the pa- 
rameter regime relevant for many-particle cooling and 
have been shown to agree well with the numerical re- 
sults. Comparison with the analogous expressions for 
free-space Doppler cooling indicate that cavity-induced 
cooling is mainly favorable for high-finesse optical cavi- 
ties fulfilling g > K. Numerical simulations also suggest 
that the cooling efficiency is increased in cavities which 
support many degenerate optical modes. 

Using the parameters of the high-Q resonators used in 
recent cavity-QED experiments one obtains spon- 
taneously scattered photon numbers per cooling time as 
low as 5-10. Replacing the heavy rubidium or cesium 
atoms by lithium, for example, therefore yields photon 
numbers below one which would allow efficient cooling 
without spontaneous emission. Similarly, only the light- 
est molecules form possible candidates for this cooling 
scheme, with the further difhculty of the reduced molec- 
ular dipole moment between any given ground and ex- 
cited state. In any case, according to the small size of 
these resonators only very few particles can be trapped 
and cooled simultaneously. 

With larger cavities the efficiency of cavity-induced 
cooling falls behind that of Doppler cooling. However, 
these systems are well suited to demonstrate the cooling 
effect on large ensembles by itself. For example, assum- 
ing sodium atoms, an atomic saturation of s = 0.1, a 
cavity decay rate k of the order of one MHz and a cou- 
pling constant g of about 100 kHz [Q, cooling times of 
the order of a ms should be achievable. 

An alternative way to reach small mode volumes com- 
bined with large finesse is the use of evanescent wave 
fields generated by total internal reflection inside mi- 
crooptical devices. Examples range from optical micro- 
spheres, thin high index surface layers or optical band 
gap guides on surfaces (integrated optics). Such devices 
look rather promising for the realization of dissipative 
walls or cooling schemes in microtraps. 



In this paper we have investigated the possibility of 
cooling a cloud of neutral particles by a cavity-enhanced 
cooling scheme. Numerical simulations of the stochastic 
differential equations describing the situation of many 
two-level atoms confined in a single-mode standing-wave 
optical resonator show that the friction force and there- 
fore the cooling time is of the same order of magnitude as 
for a single atom for the same system parameters. How- 
ever, the operating regime which is ideal for one atom 
cannot be used for many atoms, so that the best achiev- 
able cooling time increases approximately linearly with 



'^Note that for consistency with the notations used here, the 
mode functions and the pump strengths used in Ref. have 
to be multiplied by a factor of l/V^. 
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APPENDIX A: CORRELATED MOMENTUM 
AND CAVITY FLUCTUATIONS 

In the following we give the full expressions for the 
correlated noise terms introduced in Sec. |l|. The fluctu- 
ations induced by spontaneous emissions read 

^pspont ^ j^^ ^2j\£{Xn)\^ dW„ , (Al) 
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(A2) 

(A3) 



where fco = LOp/c. The noise terms according to rcscat- 
tering of photons within the cavity modes are given by 



7 Ar.ind 
1 i.ind 



dA 



k.n 



+ ^lm^^Vne'^-'^^ dW-, (A4) 



where u„ and 0„ are defined by 

V„/ (.'/;„) 
-j/fc(a;„)/2 

fk{Xn)/2 



Vnf * (s^n)^ 



|V„£:(a;„)|2- 



(A5) 



Here dH^„ is a three dimensional, all other dW^ are one- 
dimensional real Gaussian stochastic variables of mean 
zero and variance one. 
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